Magic angle effects in the interlayer magnetoresistance of quasi-one-dimensional 

metals due to interchain incoherence 
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The dependence of the magnetoresistance of quasi-one-dimensional metals on the direction of the 
magnetic field show dips when the field is tilted at the so called magic angles determined by the 
structural dimensions of the materials. There is currently no accepted explanation for these magic 
angle effects. We present a possible explanation. Our model is based on the assumption that, the 
intralayer transport in the second most conducting direction has a small contribution from incoherent 
electrons. This incoherence is modelled by a small uncertainty in momentum perpendicular to the 
most conducting (chain) direction. Our model predicts the magic angles seen in interlayer transport 
measurements for different orientations of the field. We compare our results to predictions by other 
models and to experiment. 

PACS numbers: 72.15.Gd, 74.70.Kn, 72.10.Bg, 73.90.+f 
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I. INTRODUCTION 

There is a fundamental relation between quantum co- 
herence of electronic properties and transport proper- 
ties in strongly correlated metalsi&2£&. Scattering of 
electrons affects the transport, but also blurs informa- 
tion about the momentum of the electron, and therefore 
changes the coherence of the electrons, or quasi-particle 
excitations. Generally, the effect of strong electronic cor- 
relations and incoherent excitations are enhanced in sys- 
tems of reduced dimensionality. A striking example of 
this are Luttinger liquids in one dimension. The quasi 
one-dimensional Bechgaard salts (TMTST) 2 X (X=PF 6 , 
CIO4, NO3, . . .) show a rich phase diagram ranging from 
field induced spin density waves to insulators and super- 
conductors, depending on pressure and anion 6,7 . The 
structures are highly anisotropic, and show interesting 
features as a magnetic field is applied. 

Lebed& predicted that resistance maxima would oc- 
cur when orbits along directions the crystal are com- 
mensurate with the applied field at the so called "magic 
angles" (MA) where tan# = lb/c, where 9 is the an- 
gle between the magnetic filed, tilted in the (y, z)-plane, 
and the least conducting direction, z, b and c are lat- 
tice constants, and I is an integer. The MA were later 
discovered^, but not as maxima but as dips in the an- 
gular dependence of the magnetoresistance. MA effects 
are also seer*. 10 in the (DMET-TSeF) 2 X family where 
X=AuCl2, Aula. The theory was later modified to ex- 
plain why dips should be foundii. The idea presented 
is that periodic motion is induced at the MA and, pro- 
vided there is even a small overlap in the direction of the 
applied field, the electron-electron interaction becomes 
more two dimensional which would produce a dip in the 
magnetoresistance (MR) at the MA. Alternative ideas 
and explanations have since then appeared in the litera- 
turei 12 ! 13 ! 14 ! 15 ! 16 : 17 ! 18 Most of them are based on a small 
overlap of electronic wave- functions in the direction of 
the magnetic field, but suggestions based on a Luttinger 
liquid approach have also been mad o 19 i 20 . The theory 



developed by Osada, Kagoshima and Miuraia captures 
many details of the experimental data 21 . However, that 
theory requires the existence of very long-range hopping 
(for example, the second nearest neighbor hopping inte- 
gral in a tight-binding model is of the same order as the 
next-nearest neighbor integral). Further, if B x = 0, x 
being the most conducting direction, i.e., along the one- 
dimensional chain of molecules, the theory predicts that 
there would be no MA seen in the interlayer conductiv- 
ity o zz . Further, it does not explain the dip in the MR 
when the magnetic field lies in the plane (9 = 90°). The 
data for in-plane magnetic field is affected by the fact 
that the sample is superconducting and the upper criti- 
cal field H C2 for an in-plane magnetic field is quite large^. 
The consensus is that there is no accepted theory behind 
the appearance of the MA. The experimental situation is 
also unclear at the moment. Some groups report that the 
MR has the same behavior in all directions of the cur- 
ren ^ j 2i i 23 4 24^ Whereas other experiments disagree 2 ^, and 
claim that, due to defects in the crystals, the electrons 
are forced to travel in one direction past the defects, so 
that there can be a contribution from, e.g., the resistiv- 
ity in the ^-direction to the resistivity in the z-direction, 
producing an apparently similar angular behaviour of p xx 
and p zz . 

The MA effects are also seen in torque measurements, 
as measured by Naughton et al~. This has been dis- 
cussed theoretically by Yakovenko 26 . Since the torque 
can be related to the free energy, it is likely that MA ef- 
fects reflects the ground state electronic properties of the 
material. Further, a big Nernst signal has been detected 
at the MA2£, and has been discussed theoretical!}^. 

The crystal is oriented so that x is the most conducting 
direction, followed by the y-direction. The (x, y)-plane 
defines the layered structure. Typically, the hopping in 
the three different directions are estimated to be of the 
orderiS t x :t y :t z ~ 2000^ : 200K : 10K. 

Here we present an alternative explanation of the oc- 
currence of the MA. Our physical picture is the following. 
The strongly anisotropic structure of the material affects 
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the coherence of the particles in the crystal, as well as 
increasing the effect of the electron correlations^S^i. In 
a previous paper, we discussed a model for transport 
in layered materials based on a coherence-incoherence 
crossover as a function of temperature 4 . Along the re- 
direction the motion is assumed to be coherent. The least 
conducting z-direction is assumed to be incoherent, and 
in the y-direction, the motion is predominantly coherent, 
but with a small incoherent contribution. We will show 
that the loss of coherence in the y-direction is directly re- 
sponsible for the MA in the conductance measured along 
the z-direction. Even a small amount of incoherence gives 
rise to a sizable effect seen at the MA. 



function we get that the conductivity is given by22 

Ozz — 

+G 1 -{r',r,E F )G 2+ {r,r',E F )] , 

(2) 

where G 1+ (r, r', E F ) denotes the electronic Green func- 
tion (GF) within a single layer. Here, L x L y are the di- 
mensions of the sample in the x-, and y-direction respec- 
tively. There is an indirect dependence on the distance 
between the layers in t±. We stress that this is a very 
general expression and contains all the many-body effects 
within each layer. 



II. MODEL 



We modell the system as a quasi-one-dimensional 
metal. We introduce coordinates such that a is the lat- 
tice spacing in the x-direction, b in the y, and c in the 
z, the layers lie in the (x, y)-plane. Due to the layered 
structure, the Hamiltonian is divided into intralayer and 
interlayer contributions 



H = Tin + Tt±, 



(1) 



where Tin describes the 2D (x, y)-layer and includes all 
many-body interactions within each layer, and H.± = 
t± X)<i j>{ c l c j + h.c.) describes the tunneling between 
nearest neighbors in the z-direction. Because of the lay- 
ered crystal structure, we assume that Coulomb corre- 
lations between the layers are small, and the separation 
is valid. Later, we will further specify for quasi-one- 
dimensional systems. If we have a magnetic field in the 
(y, z)-plane the vector potential, A, for the magnetic field 
B = (0,B y ,B z ) = (0, B sin 9, B cos 9). In the Landau 
gauge A is 



A=(zB y ,xB z ,0). 

We are going to study transport in the z-direction, i.e., 
transport between the anisotropic two-dimensional lay- 
ers. Let us consider two adjacent layers. The vector 
potential in the two layers are not equal but differ by 
a gauge transformation A2 = A\ + VA, 1 and 2 indi- 
cate the layer, and A = cB y x. At small bias we can use 
linear response theory to calculate the current between 
the layers. At low temperatures, only electrons at the 
Fermi-energy contribute to the conductivity in the least 
conducting direction, and it can be written as a function 
of only the in-plane Green functions^ due to the sep- 
aration of intralayer and interlayer contributions in the 
Hamiltonian. Separating the current-current correlation 



A. Non-interacting Green function for 
quasi-one-dimensional materials in a magnetic field 

Let us now look at the Hamiltonian in the absence of 
electron-electron interactions. We assume that the spec- 
tra in the most conducting direction can be linearized. 
Then, the Hamiltonian for a layer in a tilted magnetic 
field is (sec Rcf. 32) 

Wj| = av F {— ihd x +ezB sin 9)— 2t y cos [b(—ihd y — exB cos9)} 

(3) 

where a = ±1 denotes which sheet of the Fermi surface 
the electron is on. v F is the Fermi velocity, and t y is the 
interchain hopping-integral. The wavefunction is written 
as 



where 



-et 



k x x + k y y — aX sin(k y b — qx) 

(4) 



ebB cos 9 



q = 



ujb_ 

v F 



U!q cos 9 
v F 



(5) 



ojb is the frequency at which electrons traverse the quasi- 
one-dimensional sheets of the Fermi surface^ and 



A 



2f 1, 



ebvpB cost 



(6) 



is the is the wavelength of the real space oscillations of the 
electron trajectories on the Fermi surface^. In a magnetic 
field the electron dispersion relation is independent of t y , 



e a (k x ,k y ) = ahk x v F . 



(7) 



All energies are relative to the Fermi energy. The GF 
can be calculated in a way similar to the one in Ref. 13^, 
to give: 



Gl + (r,r',E) 



■— V Qe , l fc »(!'-» , )+ aJl ]e 1 ^r (l;+,r) 



(8) 
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where 



C. Interlayer conductivity 



L = sin(AyJ) — qx ) — sin(Ay5 — qx), 

and r is the electron scattering rate. The GF for the 
second layer differs by a gauge factor, e^ A ( r ~ r ' and is 



All 



(B+eu F cB sin 0+ir) 



(9) 



B. Green function containing incoherence 

We now allow for the possibility that the motion in 
the interchain direction can be incoherent. The incoher- 
ence might come from polaron formation, strong electron- 
electron correlation, or any other many-body effect. In 
a formulation in terms of GFs a possible ansatz for the 
effect of incoherence is that the non-interacting GF is 
multiplied by a y-dependent factor 



G(r,r',T)=G (r,r',T)a(y-y'), 



(10) 



where a(y — y') depends on the process by which coher- 
ence is lost. The validity of this special form of GF can be 
seen for polarons in ,e.g., Refs. 1 3331341 and for electron- 
electron interaction, in ,e.g., Ref.Ui In a 2D strongly cor- 
related model using the slave-boson approach 1 , the elec- 
tronic GF factorizes into G — GbGf, where Gf is the 
free Fermion GF, and Gg(r, r') = exp (— |r — r'| 2 /m^T) 
is the bosonic GF containing the correlations, tub is the 
mass of the accompanying boson and T the tempera- 
ture. If the 2D-lattice is anisotropic (i.e., weakly coupled 
chains) the effect from the bosonic part will be even more 
pronounced. In a previous paper we studied transport in 
layered materials of polarons^. For this case the GF con- 
tains two parts, one coherent, describing band motion of 
electrons weakly scattered by the phonons, and one inco- 
herent, where localized polarons hop between sites. For 
the case of polarons Eq. (|10fl is valid 4 ' —. Here, we do not 
specify the process responsible for the loss of coherence, 
but will just assume the general form given in Eq. (|10fl . 
The process involved in Eq. (|l(Jfl is the following. When 
the electron moves in the (x, y)-layer the k x momentum is 
conserved. Hence, there is no x-dependence in the term 
describing the incoherent contribution, o{y — y'). Instead 
it describes the change in momentum in the y-direction 
as the particle jumps between y and y' . The change in 
momentum is Sk y , which will be centered around zero so 
that for most of the time k y is unchanged. If the proposed 
form for the GF is correct, it could be visible in angle re- 
solved photoemission spectra, which measures the spec- 
tral densitji&SiSi. Later we will demonstrate that even 
a very small incoherent term gives rise to observable MA 
effects. 



Using the GFs, Eq. © and Eq. ©, in Eq. ©, and the 
incoherence factor, a(y — y'), we get a general expression 
for the conductivity: 



[dr [ dr' Y e 4^-V 
hirL x L y J J _ f-r 



a,k y ,k yl 

xkfe-y')|V 2 ^F r (e lS 



-iS\ 



)(y- y ')+ a \(L+L')] 



(11) 



where S — ecB j> m 6 | x — x'\ is the change in gauge poten- 
tial associated with interlayer transport. The summation 
over a, the two Fermi sheets the electrons moves on, can 
be done and simplified This can be simplified to 



\(L+L' 



2cos{4Asin[(A a - Ay)6/2] cos[{k y + Ay)6/2]} 

x cos {4Asin[(fc y - Ay) 6/2] cos{(k y + Ay)6/2 - q/2(x - x')]} 

+2sin{4Asin[(fc y - Ay)6/2] cos[(A y + Ay)6/2]} 

x sin{4Asin[(fc y - Ay)6/2] cos[(A y + Ay)6/2 - q/2(x - x')]} . 



Here, we introduce new variables, k y 



k y i = k + , x — x' = x + x' = x + , y — y' = y_, and 
y + y' = y+ . We can then perform the integral over x+ 
to give L x , and the integral over y + to give L y . We now 
use the representation of the trigonometric functions in 
terms of Bessel functions 

cosL4cos(A+6/2 - A)] 

J [A] + 2y(-l) k J 2k [A] cos[2A(A+6/2 - A)], 

k—l 

sin[Acos(fc+6/2 - A)] = 

OO 

2y{~\) k J 2k+1 [A] cos[(2A + l)(M/2 - A )L 

k=0 

where Ji is a Bessel function of order /. The summa- 
tion over fc_|_ can now be done by transforming it into an 
integral and we get 



Ae 2 t]_c 



Kb 



1=0 



dfc_ 



dx^ 



nv F cos 



/ ecB sm9 



4Asin 



cos(Z<ja;_)/(fc_), 



(12) 



where we introduced the distribution function 
/(*_)= / dy_e«- fe -|a(y_)| 2 , 



(13) 



describing the spread (incoherence) in the (interchain) y- 
direction. The final step is the integration in x_, which 
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gives us the final expression 

oo 

a- gz (o) = <jq 



l=—oo 



r 2 + e 2 v F B 2 (bl cos (9 -c sin (2) 2 







/ dk-Ji 


["-(¥)] 





/(*-), (14) 



where we defined the conductivity in zero field, <jq = 

^IrpbT • CSt is the main result of this paper. This 
expression can be directly compared with those derived 
by other authors for alternate theoriesi^i^. The MA ap- 
pears as peaks in a zz (dips in the MR), when the de- 
nominator has a minima. This will occur at angles when 



tan 6 = -I, 

c 



(15) 



i.e., at the MA. 

Recall that the function indicates the amount of 

incoherence in the y-direction. If we have coherent par- 
ticles in the y-direction, then, k y is always conserved so 
that ky> — k y , and the distribution will be a delta func- 
tion f(k-) = 5(k-). The sum over the Bessel functions 
collapses to only the I = term, and the result is 



a zz (B, 6) = (T 



r 2 + (ecv F B sin I 



\2 ' 



(16) 



This agrees with the result from regular Boltzmann 
transport theory^, and the MA effects are not seen. 

If an incoherent term is present we will have some 
spread in fc_. To illustrate this we use f(k-) — 



V2TTk 

the y-direction follow 



meaning that the averaged momentum in 



((ky ky/) 2 } — 



(17) 



f(k—) has the property that it becomes a delta function 
if fco — > 0, i.e., when the quasi-particles in the y-direction 
are coherent. The momentum in the x-direction is con- 
served, k x i = k x . We stress that the effects we are dis- 
cussing are not sensitive to the particular form of f(k—) 
used, since it is an integrated quantity, kob is a measure 
of how poorly the quasiparticle wave-vector is defined 
in the interchain direction. The electrons are coherent 
in the y-direction of the order of k$ , meaning that if, 
say, k$b = 0.01, then the electrons are coherent on the 
order of 100 lattice constants in the y-direction. Thus, 
a value used below kob = 0.01 still represents very well 
defined quasiparticles. A typical curve for the angular 
dependence of the interalyer magnetoresistance is shown 
in Fig. ^ The value of the other parameters, ^ and ^ 

ty 1 

are taken from typical experimental values. The decay, 
r = H/t, comes from two experiments where the scatter- 
ing time has been measured by magnetoresistance mea- 
surements and is r = 4.3ps in Ref.Hl [(TMTSF) 2 C10 4 at 
T=0.5K and ambient pressure] giving T — 0.15meV, 



and r = 6.3ps in Ref. Hi [(TMTSF) 2 PF 6 at T=0.32K 
and 8.2 kbar] giving V = O.lOmeV. The magnetic fre- 
quency, wo = ebvpB j is given by the Fermi velocity, 
vp = 0.2Mm/s in Ref. ll2t and is equal to 1.08meV when 
the magnetic field 7T and 6=7.711 A2£. The hopping pa- 
rameter in the y-direction, t y is give n as 31meV in Ref.l39l 
and Ref. IT2T but 12meV in Ref. I38L In our numerical ex- 
amples we use: ^ = 10, y 2 - =0.1. The results are not 

1 ty 

that strongly dependent on the choice of these values, 
only the amplitude of the MA dips change. Here we have 
to point out that according to the experiments^ there 
should be a dip when 8 = 90°, which is absent in our 
theory (see Fig. QJ. This dip occurs when B is parallel 
to the layers, and is therefore not a MA, and can not 
be described by our theory. As described in the intro- 
duction, it may be connected with the proximity to the 
superconducting state for the in-plane magnetic field^. 




40 

(degree) 

FIG. 1: Inter layer magnetoresistance as a function of tilt 
angle, 8, in the y — z-plane. Even for a very small incoher- 
ent hopping between the chains of molecules the magic angle 
effect are clearly seen. The parameter kob is a dimensionless 
parameter describing spread in the distribution of momentum 
as the particle tunnels between the chains, kob — means full 
coherence, i.e., a delta-function distribution of fc H -values. The 
magnetic frequency, u>o = ebvpB, is the frequency at which 
the electrons traverse the open sheets of the Fermi surface. p° 
is the resistivity in zero field, we used 6 = c and ljo/F = 10. 
We also included, as a comparison, the result when no inco- 
herence is present (kob = 0), given by Eq. Ilbt in the text. 

Note that by comparing our theory to the one by Os- 
adai^ (which assumes non-interacting electrons) the in- 
coherent term in the y-direction has a similar effect as a 
magnetic field in the x-direction. In particular we have 



B , 



hko 



(18) 



giving B x ~ 6.3T if we use fco^ — 0.01, . We see that even 
a very small incoherent part, k$b = 0.01, corresponds 
to a relatively large fluctuating field in the ir-direction 
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B x ~ 6.3T. Thus, the larger the incoherence is (larger 
kob) the larger the corresponding effective field in the re- 
direction is, and the larger the MA dips in the MR are. 
This is consistent with the experimental result by Lee 
and NaughtonSi, where an increasing ^-component of the 
magnetic field increased the size of the MR oscillations 
at the MA. 



III. B IN THE (x, z)-PLANE 

If we instead apply the magnetic field in the (x, z)- 
plane the vector potential will be 

A= {0,xB z - zB x ,0). 

The derivation is very similar to the one presented above 
with the only difference that the gauge potential does 
not have any component depending on \x — x'\, but now 
depends on y — y' instead. The result is that the inte- 
gral over X- is simpler, but the integral over ?/_ has an 
additional factor. This factor can be absorbed in the y_ 
integral, the final result is 



(TO 



;=-c 







fdk-Ji 


M¥)] 





r 2 



T 2 + (ev F blB z ) 



(19) 



where 

g(k-) = J dy_e^-( fc 
The parameter A is 



Hv-)? = f [k 



21 „ 



ecB x 

h 

(20) 

The so called Danner-Kang- 



ebvpB 

Chaikin oscillations^ are observed provided that 



ecB r 



> k 



(21) 
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FIG. 2: Interlayer magnetoresistance as a function of the 
magnetic field direction in the x — z-plane. 9 is the angle 
between the field and the z-axis. The upper panel shows a 
numerical calculation of the so called Danner-Kang-Chaikin 
oscillations^, from Eq. I19H in the text. The theoretical curve 
can be compared with Fig. 1 from Ref. |38| shown in the lower 
panel, with experiments done on (TMTSF^CICU at ambient 
pressure and T=0.5K. The dip around zero degree below 3T 
is due to the sample becoming superconducting. Note that 8 
denotes the angle between the magnetic field and the i-axis. 



We used b - 
at B = IT. 



c, kob = 0.001, with ujo/T = 10 and uo/ty 



0.1 



where kgb is the incoherence parameter. In Fig. [2] we 
compare the resulting resistivity (l/er zz ) from Eq. i|19|) 
with an experimental curved. We did not adapt the pa- 
rameters to the experiment, but just want to illustrate 
that this type of oscillations do appear in the theory 
presented. Note that we have used a smaller value for 
the incoherence parameter kob — 0.001, compared to the 
value used in Fig.^ This is justified by the fact that the 
experiment we compare with is performed for the CIO4 
compound and the oscillations in the y — z-plane are not 
as visible^ as for the PFg compound indicating a smaller 
incoherence factor. 



IV. B IN THE (x, y, z)-PLANE 

Combining the results from the calculations above we 
can get an expression for a field in a general direction, 



(B x ,B y ,B z ). We get 

OO 

Vzz{Q,<t>) = ^0 ^2 

I — — OO 

dk-Ji 



r 2 + e 2 Vp(blB z - cB y ) 2 
1 ,<>■ ii \ 

4Asin 



fc_6 
2 



f[k-- 



ecB 7 



(22) 



note that A is a function of B z . In Fig. [3] we compare 
results from this expression with the experimental results 
of Lee and NaughtonS, by identifying the angles defined 
in Fig. |3 as follows, 



B x = B cos 6 cos < 
By = B cos 9 sin <j 
B 7 = B sin 9 



(23) 



6 



where the definition of and <j) follows Ref. |3^, (see the 
upper panel in Fig.|3J). As the angle 9 between the (x, y)- 
plane and the direction of the field is increased, the os- 
cillations start to appear. The similarities to Fig. 4 in 
Ref. |3^ are striking. Changing the parameters in the 
model does not change the general features of this plot. 

V. DISCUSSION 

In summary we have presented an explanation in terms 
of many-body effects of the appearance of magic angle 
effects in the interlayer magnetoresistance. The MA ap- 
pears naturally from, even a small, incoherent contribu- 
tion to the inter-chain hopping. The hopping in the most 
conducting direction is assumed to be coherent, and in 
the least conducting direction incoherent. Momentum 
can change in the direction between the one-dimensional 
chain of molecules. This is described by a distribu- 
tion function which is centered around zero, letting most 



quasi-particles retain their momentum when hopping. 
We used an explicit form of the interlayer Green function, 
which can be directly observed in a angle resolved photoe- 
mission spectra. Unlike present explanationsiSiSii^l, the 
theory does not assume any long distance hopping be- 
tween non- adjacent quasi-one-dimensional molecules in 
different layers, where the overlap is quite small, only 
a nearest neighbor interlayer overlap. The shape of the 
Fermi surface is not affected by the incoherence. Numer- 
ical calculations produce results similar to experimental 
results. 
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FIG. 3: Interlayer magnetoresistivity versus j/— z-plane angle, 
a, denned via tana = sin</>/ tan# (see top figure). The mid- 
dle panel shows the result from our numerical calculation of 
conductivity using Eg. 1221 Modulations appear at the magic 
angles as the angle a is increased. We used b — 7.58lA and 
c = 13.264A. The other parameters used are kob — 0.1, 
tJo/T — 10 and uJo/t y — 0.1. The theoretical curve can be 
compared with Fig. 4 from Ref. shown in the lower panel. 
This is an experiment done at 0.32K on (TMTSF) 2 PF 6 with 



